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ABSTRACT 

We calculate orbits for the Milky Way dwarf galaxies with proper motions, and com- 
pare these to subhalo orbits in a high resolution cosmological simulation. We use 
the simulation data to assess how well orbits may be recovered in the face of mea- 
surement errors, a time varying triaxial gravitational potential, and satellite-satellite 
interactions. For present measurement uncertainties, we recover the apocentre r^, and 
pericentre to ~ 40%. With improved data from the Gaia satellite we should be able 
to recover ra and rp to ~ 14%, respectively. However, recovering the 3D positions and 
orbital phase of satellites over several orbits is more challenging. This owes primarily 
to the non-sphericity of the potential and satellite interactions during group infall. 
Dynamical friction, satellite mass loss and the mass evolution of the main halo play a 
more minor role in the uncertainties. 

We apply our technique to nine Milky Way dwarfs with observed proper motions. 
We show that their mean apocentre is lower than the mean of the most massive 
subhalos in our cosmological simulation, but consistent with the most massive subhalos 
that form before z = 10. This lends further support to the idea that the Milky Way's 
dwarfs formed before reionisation. 
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O 1 INTRODUCTION 



The dwarf galaxies of the Milky Way provide a unique win- 
dow onto galaxy formation on the smallest scales. Their 
close proximity allows us to resolve individual stars, obtain- 
ing detailed star formation hist ories, mass m easurements, 
proper motions, and distances (jMateol [l998l ). Increasingly 
accurate proper motions present th e possibility of reliable 



orbit calculations for the dwarf s fe.g.lSchweitzer et al.l 
Schweitzer fc CudworthI 1 19961: iPiatek et all |2002|. 



1995 



200a 



Dinescu et al.1 |2004 iPiatek et al.l l2005l : iKallivavalil et all 
20061 : ICosta et al.ll2009l ). We may then explore links between 



environment and f ormation history, e.g . the morphology- 
distance relation [(Ivan den Berghl 1 19941 ) although recent 
discovery of three dSph at large distances from An- 
dromeda and the Milky Way calls this into question (e.g. 
iMcConnachie et al.l |2008| )] . With accurate 3D determina- 
tions of their orbits, we may then also determine if the 
dwar f s fell into the Milky Way in groups (jP'Onghia fc Lakj 
l2008l : iLi fc Helmil l2008l: iKUment owski et al.' ' 20091) or su f- 
fered past interactions (|Sales et a l. 2007; Matco ~et al.ll200§ ). 

In this paper, we derive orbits for the Milky Way's 
dwarfs with proper motion data, by integrating their mo- 
tion backwards in a static gravitational potential. We es- 
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timate the effects of several systematic errors by using a 
similar integration of subhalos in the high resolution Via 
Lactea I cosmological simula tion of a Milky Way analogue 
(VLl; iDiemand etallbOOTah . We assess the impact of mea- 
surement errors by adding realistic errors to the subhalos' 
redshift z — phase space position, and tracing their orbits 
back in a fixed potential. Our fiducial model integrates or- 
bits in a fixed spherical potential; however we also examine 
the effects of a triaxial potential, mass evolution of the main 
halo, dynamical friction, and mass loss from the satellites. 
We devise three measures of quality for tracing back orbits, 
of increasing difficultly. These are the comparison of true 
and recovered: 

(i) last pericentre Vp and apocentre Va 

(ii) Tp, ra and the orbital period t backwards in time over 
A'^ orbits. 

(iii) 3D pericentre rp and apocentre ra backwards in time 
over A'^ orbits. 

The first of these three tells us the current orbit of the 
dwarf for comparison with cosmological simulations. It rep- 
resents the minimum useful information we might obtain 
from a dwarf orbit. The second tells us when the dwarf en- 
countered the Milky Way, and how close it came. This allows 
us to compare orbits with star formation histories. The third 
tells us the full 3D position of the dwarf as a function of time. 
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This tells us if dwarfs ever interacted via group infall or fast 
flyby. 

This paper is organised as follows. In !j2] we exam- 
ine whether orbits can be recovered at levels (i), (ii) and 
(iii) within the model limitations and what data quality 
is required for this. In ^ we apply our method to nine 
Milky Way dwarf galaxies with measured radial velocities 
and proper motions. Finally, in U we present our conclu- 
sions. 



2 TESTING THE METHOD 

In this section, we use the Via Lactea I (VLl) sim ulation of 
a Milky Way mass galaxy (jPiemand et al.' 2007b'), to deter- 
mine how well we recover satellite orbits in the face of mea- 
surement errors and a time varying gravitational potential. 
We extract from the simulation three sets of subhalo^: the 
50 most massive today (zo°)i the 50 most massive before red- 
shift z — 10.59 (^lo) and the 50 most massive before redshift 
z = 10.59, taking depletion by a disc into account (zfo (rd) 
In all cases, we include only subhalos with mass 
and distance to the centre of the main halo r < 150 kpc at 
redshift z = 0. This represents the mass and radius range 
where we find the Milky Way dwarfs (c.f. Table [T] in sec- 
tion i ]3.ip . In the disc depleted sample we exclude all orbits 



having pericentres < before extracting the sample, mo- 
tivated by recent work bv lD'Onghia et al] (|2009l '). We use 50 
subhalos to create a sensible sized sample that likely corre- 
sponds to the expected number that might be observed with 
full sky coverage at a depth comparable to the Sloan Dig- 
ital Sky Survey JSimon fc Gehall2007l : IKoposov et ai]|20o4 
IWalsh et al.ll2009l v 

We take the present day phase space position of these 
satellites and integrate them backwards in time to com- 
pare our derived orbits with the true VLl orbits. We use 
a Leapfrog time integrator with adaptive time stepping to 
adjust for the hig her resolution required at pericentre (e.g. 
iPress et al.lll993 ). We scale each time step b y a fraction 
of the instantaneous orbital time, similarly to IZemp et al.l 
(|2007f ). The apocentre and pericentre, rp of each orbit 



(and therefore the period t) are then recovered by search- 
ing for a sign change in dr/dt, where r is the distance from 
the satellite to its host galaxy. For the VLl data, we use a 
more sophisticated algorithm, that only characterises global 
extrema in the orbit around the main halo and ignores lo- 
cal extrema caused by satellite-satellite interactions or the 
satellite behaviour before falling into the main halo. We ex- 
plicitly exclude the cosmological turn around as part of the 
real orbit. Note that for both the simulation as well as the 
recovered orbits, we determine half periods ti/2, i.e. the tem- 
poral distance between two orbital extrema. Hence, we only 
include orbits with 2 orbital extrema. By this criteria any 
satellites with long periods are excluded, which are falling 
in for the first time. Whenever we write period instead of 
half period, we mean the half period multiplied by two. 
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Figure 1. Evolution of the mass and scale radius of the VLl main 
halo (solid line). The scale radius is more difficult to constrain for 
high redshifts. The dashed line is the best fitting power law going 
back 9 Gyrs in time. 



2.1 The orbit recovery models 

We consider the following models for our orbit integration: 

• The fiducial model (F) uses a static, spherical, NFW 
potential (|Navarro et al.|[l996l ). In F, we take parameters 
from the best fit to the VLl main halo at redshift 2 = 0: 
M200 = 1.77 X 10^^ M q, R200 = 389 kpc and Vc = 24.6 kpc 
(jPiemand et al.ll2007al ). 

• The dynamical friction model (DF) uses the same 
potential as in model F with dynamical friction forces added 
(|Chandrasekhailll943f). We use a radially varying Coulomb 
logarithm as in lHashimoto et al. We also explore the 
impact of mass loss of the infalling satellites, and the mass 
growth of the VLl halo with time (see Figure [TJ. We fit 
simple functions to the subhalo mass l oss his t ory, a nd VLl 
main halo mass growth with time as in IZhad (|2004l l. 

• The tri axial model ( T): T his u ses a triaxial NFW po - 
tential as in lKuhlen et all (|2007l ') and ICuedes et all |200i). 
We use axis ratios: q = 0.83 and s = 0.8, as measured for 
the VLl main halo at outer radii at z = 0. 

• The double halo mass model (2Af): This uses the 
potential as in model F, but doubling the mass of the main 
halo (doubling the scale radius has a negligible effect). 



These different models allow us to assess the importance 
of systematic effects on our orbit recovery for real Milky Way 
dwarf data. In practice, only model F can be realistically 
applied to real data since we do not know the mass loss his- 
tory of the Milky Way dwarfs, nor the mass growth history 
or shape of the Milky Way halo. We use models DF, T and 
2M to explore the systematic impact of this poor knowledge 
on our orbit recovery. 



2.2 Results 

In this section, we measure how well we can recover the 
orbits of subhalos in our VLl data set. We use each of the 
the recovery models presented in H2.ll with and without 
measurement errors and assess our results using criteria (i), 
(ii) and (iii) in 



2.2.1 Perfect data 



^ We refer to subhalos as 'satellites' interchangeably. 

^ See the data at |http://www.ucolick.org/~diemand/vl/] 



To assess the model systematics, we first consider the case 
without position and velocity measurement uncertainties. 
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Figure 2. Orbit recovery in models F, DF, T and 2M with no measurement errors. The coloured bars show the percentage of subhalos 
with fractional error Q = — 0.1 (blue), 0.1 — 0.3 (cyan), 0.3 — 0.5 (green), 0.5 — 0.8 (red) and > 0.8 (magenta) over N orbits backwards 
in time. In our Fiducial model (_F), 20% of satellites have their most recent pericentre recovered to better than 10%, while 56% have 
their most recent apocentre recovered to better than 10%. 
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Figure 3. As Figure[2l but for the 3D pericentre (left) and 3D apocentre (right). 



Figures [2] [3] and |4] show our recovery for each of the models 
F, DF, T and 2M for the 'ID' pericentres rp and apocentres 
ra, the 3D pericentres rp and apocentres ra and the half 
period ti/2 over A'^ orbits. The coloured bars show the per- 
centage of subhalos with fractional error Q — — 0.1 (blue), 
0.1 - 0.3 (cyan), 0.3 - 0.5 (green), 0.5 - 0.8 (red) and > 0.8 
(magenta). We define the relative error Q as: 

Q^^r_p-j^ . Q^(30) ^ ^.ZlP:^ (1) 

for the 'ID' pericentre (left) and 3D pericentre (right), where 
Vp.t is the true pericentre and Vp is the recovered pericentre. 
Q is defined similarly for the other orbit quantities Va, fa 
and ti/2- Note that in comparing Vp and rp,t, the limited 
time resolution of the VLl simulation can lead to an addi- 
tional error (that causes us to always over estimate rp^t)- 



We tested the effect of this, finding that only orbits with 
rp < 20kpc are affected. The additional error due to finite 
time resolution is small compared to our model uncertain- 
ties. 



2.2.2 Recovering the ID pericentre and apocentre 

From Figure [2l we see that the last apocentre is well recov- 
ered in our fiducial model F. No satellites have their last 
apocentre recovered with greater than 30% error. This is 
perhaps to be expected. Subhalos are most likely to be at 
apocentre now and so their present position already gives 
a reasonable constraint on r^. The DF and T models both 
improve the recovery for the most recent apocentre. This 
suggests that the primary reason for our error in the most 
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recent apocentre recovery is a combination of wrong halo 
shape and dynamical friction, where the halo shape is more 
significant. These effects will be swamped by other effects 
for the apocentre recovery further back in time (c.f. I2.2.3p . 
By contrast, model 2M gives a poor performance, with only 
~ 20% of the subhalos having better than 10% error on the 
last apocentre. This demonstrates that it is critical to have 
a good estimate of the mass of the Milky Way interior to 
the orbits of its dwarfs. 

The last pericentres are significantly less well recovered 
than the apocentres, with only 20% having an error better 
than 10% in our model F. Here the halo shape plays a major 
role, with model T giving a significantly better performance 
(~ 40% recovered to better than 10%); dynamical friction 
(model DF) has little impact on the recovery. To better 
understand why the pericentre recovery fails, some typical 
example orbits from each Qp bin are given in Figure (5] As 
can be seen, orbits that have a well recovered last Vp are 
short period orbits. These are well recovered in model F and 
well recovered in the triaxial model T that has the correct 
current halo shape. However, long period orbits sample the 
potential over several giga years backwards in time. Here 
even our triaxial model (fit to the VLl halo at redshift z = Q) 
fails and the rp recovery is poor (see the magenta orbit). 

2.2.3 recovenng orbits over several periods: The problem 
of group infall 

Recovering orbits over several periods backwards in time is 
more challenging. The fraction of orbits with better than 
10% error in ra in our triaxial model T drops by ~ 40% by 
Na = 2; while there are no well recovered rp in our model F 
by Np = 3. The triaxial model T gives an improvement in 
the recovered rp, but no improvement in r^. Note that there 
are some fiuctuations in this trend due to small number 
statistics in the Na ~ 3 and Na = 4 bins. So individual 
orbits, that are improved in one of the models have a higher 
impact than they would have in a larger sample. 

To understand why the errors on ra grow so rapidly 
backwards in time, we give some example orbits recovered 
at Na = 2 in Figure [B] As with the rp (Figure [5]), the longer 
period orbits are more poorly recovered. 

There is a significant class of orbits - 50% - with 
0.1 < Qa < 0.3 (cyan) that are less well recovered at Na = 2. 
We explore a typical orbit in this class in more detail in 
Figure [7] The middle panel shows our orbit recovery for the 
most massive subhalo in this group in model F (black dashed 
line), T (magenta dashed line), DF (blue dash-dotted line) 
and an orbit integrated with dynamical friction, but with- 
out mass loss of the satellite nor mass evolution of the main 
halo (green dotted line) . Notice that none of our models give 
a good fit to the orbit. Accounting for dynamical friction 
and mass loss (see right panel) as well as using the triax- 
ial potential significantly overestimates the apocentre. We 
investigated this orbit further and find this satellite to be 
part of a loose group. The left panel shows the orbit of the 
satellite (black) and all subhalos that were initially close in 
phase space to this satellite (red). This grouping was found 
by determining all sate llites closer than 4rt, w here rt is the 
tidal radius defined in iDiemand et~al] (|2007al ). for at least 
six time outputs. These 'loose groups' usually break up at 
first pericentre, which causes an energy change that results 



in lower 2nd apocentre than we would expect from any of 
our integration models. 

75% of the the satellites in the cyan bin were found to be 
falling into the galaxy as part of a group. This s uggests that 
the gr oup infall statistics recently determined bv lLi fc Helmil 
(|2008l ) may be a lower bound. We defer a detailed analysis 
of the statistics of loose group infall to future work, but note 
here that it appears to be responsible for many of our poor 
orbit determinations at A''^ ^ 2. 

2.2.4 recovering the orbital phase 

In Figure |4l we show how well we recover the half orbital 
period over Ni^2 orbits assuming perfect data. Notice that 
the period is well-recovered up to even two orbits backwards 
in time. However, in determining the phase of the orbit, such 
period errors accumulate. This makes it challenging to try 
and match pericentre passages with observed star formation 
histories for the dwarfs. By N1/2 ~ 2, our typical phase error 
is 0.6 Gyrs, at Ni^2 =3 it is already 0.8 Gyr. Recall that this 
is for perfect data that have no measurement errors. With 
current proper motion errors, recovering the orbital phase is 
simply not possible (see ^2.2.6y 

2.2.5 recovering the 3D pericentre and apocentre 

Figure O shows how well we recover the 3D pericentre and 
apocentre for perfect data. Here having an accurate halo 
shape is vital. Our results in model F, even for a single or- 
bit, are poor. Even in our triaxial model T the results, while 
dramatically improved, are not encouraging. This is most 
likely due to the ra dial and temporal variations in the triax - 
iality as found by l|Kuhlen et al.ll2007l : ICuedes et al.ll2009l '). 
~20% of the most recent 3D apocentre distances Ta and 
less than 5% of the most recent 3D pericentre distances Tp 
are recovered to better than 10% accuracy. Our results sug- 
gest that, unless the Milky Way halo is close to spherical or 
axisymmetric, recovering full 3D orbits for the Milky Way 
dwarfs will not be possible even with perfect data. This will 
make it difficult to determine if any of the dwarfs fell in to 
the Milky Way together, or had past interactions with one 
another. 

2.2.6 The effect of measurement errors 

In this section, we add the effects of measurement errors. 
We consider both current typical errors, and those expected 
in the Gala era. For current errors, we use Ar — 5kpc, 
Awjos = 30km/s and Avt = 60km/s, where Ar is the dis- 
tance error and Avios,t are the line o f sight and tangen - 
tial velocity errors, respectively (see e.g. IPiatek et al.ll2007l ). 
For Gaia errors, we use Ar = 5kpc, Avios = 5k m/s and 
Avt = lOkm/s (see e.g. I Wilkinson fc EvansI [l999h . We de- 
termine the effects of measurement errors on our derived 
orbits by building an ensemble of 200 orbits drawn from the 
above distributions around an already defective orbit, where 
the defective orbit is also drawn from the same error distri- 
bution. This represents the simplest way to mimic the real 
measurement errors of the Milky Way dwarfs. 

We show results for our Fiducial model F in Figure |8l 
We plot the true (black crosses) and recovered (blue error 
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Figure 7. Group infall in the VLl simulation. Tiie left panel shows the orbits of subhalos that fall in inside a group (red and black 
lines); the most massive member is shown in black. The middle panel shows our orbit recovery for the most massive subhalo in this group 
in model F (black dashed line), T (magenta dashed line), DF (blue dash-dotted line) and an orbit integrated with dynamical friction, 
but without mass loss of the satellite nor mass evolution of the main halo (green dotted line) . The right panel shows the mass evolution 
of the most massive subhalo in the group and our best fit curve (blue dot-dashed line) used in model DF. 
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Figure 4. As Figure[2l but for half the orbital period ti/2- 

bars) pericentres rp, apocentres Va and their ratio Vp/ra- 
The black dashed line shows the mean of the true values, 
while the grey shaded band marks one sigma scatter around 
the mean. The blue dashed line is the mean of the recovered 
orbits. 

For current measurement uncertainties, the error is 
dominated by the proper motion errors (A«t). For Gaia er- 
rors, we are instead limited by model systematics. However, 
such systematics appear to average out over the whole pop- 
ulation. Although our mean recovered rp and Va are both bi- 
ased high by the proper motion, this biasing is significantly 
reduced with Gaia quality data. 

2.2.7 Orbit recovery at levels (i), (ii) and (in) 

With current proper motion errors, we are able to recover 
orbits only at level (i) (the last pericentre rp and apocentre 
ra\ see With Gaia quality data, it will be possible to 
recover orbits at level (ii) (rp, and the orbital period t 
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Figure 5. True and recovered orbits for the first pericentre with 
different Qp values. The colours correspond to the bins in Figure 
[2] Qp = 0-0.1 (blue), 0.1-0.3 (cyan), 0.3-0.5 (green), 0.5-0.8 
(red) and > 0.8 (magenta). Solid lines are the true orbit; dashed 
lines are the recovered orbits in our fiducial model F\ and dotted 
lines are the recovered orbits in our triaxial model T. 



backwards in time over N ^ 2 orbits), though the pericen- 
tres will suffer from large systematic errors if the Milky Way 
potential is triaxial. Even with perfect data, going back fur- 
ther than ~ 2 orbital periods runs into the problem of group 
infall (see i]2.2.3[) . Full 3D orbit recovery (level (iii)) will be 
extremely challenging, recovering 3D orbit data will only be 
possible if the Milky Way potential is nearly axisymmetric 
or spherical, and if it did not change significantly over the 
past ~ 8 Gyrs. 
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Figure 8. Orbit recovery in our Fiducial model F, with measurement errors. Panels show the pericentre, apocentre and peri/apo ratio 
for the 50 most massive satelhtes at ^ = in the VLl simulation. The black crosses correspond to the VLl data, the blue error bars to 
the recovered data from the orbit integration. The black dashed horizontal lines and grey shaded area correspond to the mean and the 
one sigma variation of the true values, respectively. The blue dashed line corresponds to the mean of the recovered values. 



3 APPLICATION TO THE MILKY WAY'S 
DWARFS 

In this section, we apply our orbit recovery technique to 
nine Milky Way dwarfs with observed proper motions. We 
compare our derived orbits to those of subhalos in the VLl 
simulation. With current measurement errors, we can only 
obtain a reliable estimate of the last pericentre and apocen- 
tre (see H2.2.7|l and so we focus on these orbit diagnostics. 



3.1 The data 

Table [T] summarises the dwarf data we take from the liter- 
ature. From left to right the columns show: galactic longi- 
tude I and latitude h in degrees; the distance to the sun d in 
kpc; the radial velocity Vr in km/s; the proper motion data 
jia cos(5), fj,5 in mas/yr; the V-band magnitude Mv (mag); 
and the Mass within 600 pc in units of IO^Mq. We order 
the dwarfs by their star formation history, placing them in 
three broad categories: those with mainly early star forma- 
tion (> 9Gyrs; E); those with mainly intermediate age star 
formation (3 — 9Gyrs; I); and those with any significant re- 
cent star formation (< 3Gyrs; R). Some half-categories are 
also introduced for galaxies that have mostly early star for- 



mation, but evidence for some intermediate age also (EI). 
The orbit data recovered in this paper as well as our classi- 
fication of the star formation histories and the corresponding 
references are summarised in Tabled 



3.2 Potentials 

For the Milky Way potential, we use two different mass mod- 
els from the literature. This gives us a handle on the sys- 
tematic error arising from our potential model. For our first 
model, we use the oblate potential from Law et al ] (l2005h 
consisting of a Miyamoto-Nagai disc ('Mi yamoto fc Nagail 
[l975), Hernquist spheroid |Hernauist..l990, l . and a logarith- 
mic halo. The LOS model: 



<&disc = 



<E>: 



sphere — 



phere 



r + c 



(2) 

(3) 

(4) 
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Galaxy 


I (dcg) 


b (dcg) 


d (kpc) 


Vr (km/s) 


fia cos((5) (mas/yr) 


(mas/vr) 


A'/v fmag) 


AfeooflO^Afol 

uuu L j-'j- J 


Refs. 


UMi 


105.0 


+44.8 


66 ±3 


-248 ± 2 


-0.50 ±0.17 


0.22 ±0.16 


-8.9 




1,2,12 


Draco 


86.4 


34.7 


82 ±6 


-293 ± 2 


0.6 ±0.4 


1.1 ±0.3 


-8.8 




1,3,12 


Sextans 


243.5 


+42.3 


86 ±4 


227 ±3 


0.26 ±0.41 


0.10 ±0.44 


-9.5 


0.9lgj 


1,4,12 


Sculptor 


287.5 


-83.2 


79 ±4 


108 ±3 


0.09 ±0.13 


0.02 ±0.13 


-11.1 


9 7+0.4 
^■'-0.4 


1,5,12 


Carina 


260.1 


-22.2 


101 ±5 


224 ±3 


0.22 ±0.13 


0.24 ±0.11 


-9.3 


3.41?:^ 


1,6,12 


Fornax 


237.1 


-65.7 


138 ±8 


53 ±3 


0.476 ± 0.046 


-0.360 ±0.041 


-13.2 




1,7,12 


Sagittarius 


5.6 


-14.1 


24 ± 2 


140 ±5 


-2.35 ±0.20* 


2.07 ±0.20* 


-13.4 


97+20 
^' -27 


1,8,12 


SMC 


302.8 


-44.3 


58 


175 


1.16 ±0.18 


-1.17 ±0.18 


-17.2 


10 ±2 


1,9,10,13 


LMC 


280.5 


-32.9 


49 


324 


2.03 ±0.08 


0.44 ± 0.05 


-18.6 


14 ±3 


1,10,11,14 



Table 1. Dwarf data from the literature. From left to right, the columns show galactic longitude; galactic latitude; distance to the sun; 
radial velocity ; prop er motions; V-band ma gn itude; the dynamical m a ss within 600 pc (Mg op) ; and the dat a refer en ces. Data were taken 
from: UMatco' (1998); ^Piatck ct al." ('2005'); ^'Scholz & Irwin ('1994'); -^Walker et al.' ('2008'); ^Piatck ct al.' (2006); ^'piatck ct al.' {200A 
^004); 'Piatck ct al. (2007); * Dincscu et al. (2005); ^ Kallivayalil ct al . (20 06); Koposov ct a l. (2009) ; " Kallivayalil ct al. {20oM; 
^^Istrieari et a l. (2007); ^'^ mass derived from me an oflHatzidimitriou et al. I | |T997l '): IStanimirovic et al.H 2004'): H arris Zaritskvl l |2o6a) : 
lEvans &c Howarth. ( 2008 J ; mass derived from iKim et al. I ( Il998l) . consistent with Ivan der Marel et al., (i2002i ) ; * these values are for 
in cos(6) and fii,, respectively. 
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Figure 6. True and recovered orbits for the second apocentre 
(Na = 2) with different Qa values. The colours correspond to the 
bins in Figure[2] = - 0.1 (blue), 0.1 - 0.3 (cyan). Solid lines 
are the true orbit; dashed lines are the recovered orbits in our 
fiducial model F; and dash-dotted lines are the recovered orbits 
in our dynamical friction model DF. 

with Mdisc = 1.0 X 10" Af©; a = 6.5kpc; b = 0.26kpc; 
Msphorc = 3.4 X 10^°Mq; c = 0.7 kpc; -Uhaio = 121 km/s; and 
q = 0.9. 

For our second mode l , we u se the Truncated Flat model 
from IWilkinson fc EvansI (|l999l ). The TF model: 

= — log( 1. (5) 

with values M = 1.9 x lO^^M©; and a = 170 kpc. 

These different models have quite different asymptotic 
total masses, but the mass within the orbit of the dwarfs 
(~ 150 kpc) is in reasonable accord (~ 6 x lO^^M© in the 
LOS model, and 12 x lO^M© in the TF model). The value 
for the VLl is 10 x 10^^ and well within that range. We 
can therefore directly compare our results without the need 
to rescale the simulation. 



Galaxy 


rp (kpc) 


ra (kpc) 


T(Gyrs) 


SFH 


Rofs 


UMi (L05) 


40 ±20 


90 ±20 


1.4 ±0.4 


E 


1,2 


UMi (TF) 


30± 10 


80 ±10 


1.2 ±0.2 


E 


1,2 


Draco (L05) 


90± 10 


300 ± 100 


6±2 


E 


1,3 


Draco (TF) 


74 ±6 


250 ± 80 


4±2 


E 


1,3 


Sextans (L05) 


70 ±20 


300 ± 200 


4±4 


E 


1,4 


Sextans (TF) 


60 ±20 


200 ± 100 


4±3 


E 


1,4 


Sculptor (L05) 


60± 10 


160 ± 80 


3±1 


EI 


1 


Sculptor (TF) 


60± 10 


130 ± 60 


2±1 


EI 


1 


Carina (L05) 


60 ±30 


110 ±30 


2.0 ±0.6 




1 


Carina (TF) 


50 ±30 


110 ±30 


1.8 ±0.8 




1 


Fornax (L05) 


120 ± 20 


180 ± 50 


4±1 




1,5 


Fornax (TF) 


110 ±20 


170 ± 40 


4±1 




1,5 


Sagittarius (L05) 


12 ± 1 


50 ±7 


0.7±0.1 




1 


Sagittarius (TF) 


12 ± 1 


53 ±5 


0.56 ±0.08 




1 


SMC (L05) 


57 ±5 


200 ± 100 


4±1 


R 


1,6 


SMC (TF) 


56 ±5 


200 ± 100 


3±2 


R 


1,6 


LMC (L05) 


47 ± 1 


500 ± 100 


7±2 


R 


1,7 


LMC (TF) 


47 ± 1 


400 ± 100 


7±3 


R 


1,7 



Table 2. Derived orbits for the dwarfs, along with SFHs from the 
literature. From left to right, the columns show: Galaxy name; 
pericentre in kpc; apocentre in kpc; period in Gyrs; the Star for- 
mation history classification (see il3.1|l : and the data references. 
Each galaxy appears twice to show results from the LOS poten- 
tial and the TF potential (see il3.2|l . Notice, that even though 
the apo- and pericentre are biased high with current measure- 
ment errors, the period is biased low. This is due to the bias 
towards more circular orbits. The star formation history data 



were taken from: HOolphin et a l. (2005); ^Carrera et~all (120021); 



^Aoaricio et al.l (l200lft: ' 


Lee et al.. (.2009): ^Coleman & de Jond 


(l2008f); iNoel et al.l (12009 


): iHarris & Zaritskvl (l200gf). 



3.3 Results &: discussion 

In Figure [9] the current radial distribution of the nine MW 
dwarfs with errors are shown. Overlaid are the current ra- 
dial distributions of the two samples Zq^ and zIq. In seem- 
ing contradiction to former results in th e literature (e.g. 
iKravtsov eral] 12004 iMaccio' et al] |2009| ). both distribu- 
tions are equally consistent with the current dwarf distri- 
bution. This surprising result can be explained by a pecu- 
liar phase distribution of the orbits. In Figure 1101 the radial 
distribution of Zq° and Ziq averaged over the last 2 Gyrs 
are shown. On this time scale the orbits do not change sig- 
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Figure 9. Current radii for nine Milky Way dwarfs with observed proper motions with errors. Overlaid are the mean and standard 
deviation of the values from VLl for the 50 most massive at z = (left panel) and for the 50 most massive before z = 10 (right panel), 
respectively (grey band, dashed line). 
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Figure 10. As Figure O but for the radial distributions averaged over the last 2Gyrs. 



nificantly, but their phases are averaged out. Now we see a 
distinct difTerence between the two distributions. Therefore 
the apo-/pericenter distributions of both samples need to be 
compared with the dwarfs to differentiate between different 
dwarf formation scenarios. (Note that previous works using 
the same VLl simulation that we use here did see a differ- 
ence in the radial distribution of dwarfs selected at z = 
and those selected a.t z — 10. This is because they did not 
use the same radial and mass cuts that we employ here. For 
real data, we may be lucky and see a biasing in the radial 
distribution, but to be certain to see it, we must use the 
apo/pericentre distributions.) 

In Figure [TT] we show our derived orbits for nine Milky 
Way dwarf galaxies with observed proper motions. Analo- 
gous to the VLl satellites (c.f. i]2.2.6|) . we estimate our errors 
by building an ensemble of 1000 orbits for each dwarf drawn 



from its error distribution. The error bars show the values 
for the oblate L05 potential (blue) and for the TF model 
(red). Overlaid are the mean and standard deviation of the 
values from VLl for the 50 most massive subhalos at z = 
and for the 50 most massive before z = 10, respectively 
(black dashed lines and grey bands). The blue dashed lines 
denote the recovered mean of the Zq^ and zfo distributions, 
respectively (c.f. Figure [8]). Note that these lines indicate a 
tendency, as they only represent one possible realisation of 
the subhalo distribution. The galaxies are ordered by their 
observed star formation histories as described in 

We do not model the LMC and SMC together as is 
often done in the literature (e.g. iBesla et al. I l2007l : iBekkil 
|2008| ). but investigate instead their orbits independently of 
each other. As shown in i]2.2.5l our model systematics are 
too large to determine whether or not it is better to include 
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Figure 11. Recovered orbits for nine Milky Way dwarfs with observed proper motions. The error bars show the values for the oblate 
L05 potential (blue) and for the TF model (red). Overlaid are the mean and standard deviation of the values from VLl for the 50 most 
massive at z = and for the 50 most massive before z = 10, respectively (grey band, dashed line). We also overlay the recovered mean 
as in Figure |8] (blue dashed line). Note that these lines indicate a tendency, as they only represent one possible realisation of the satellite 
distribution. The galaxies are ordered by t heir observed sta r formation histories as described in i|3.1l The black diamonds denote the 
values derived from the Sagittarius stream (|Law et al.ll2005h . 



LMC-SMC interactions in our models. Notice, that we only 
find bound orbits for the LMC because of the large Milky 
Way masses in our models. Recen t results for LMC orbits 
in a lower MW mass p otential bv iKallivavalil et al.l ()2006l ) 
and iBesla et all (|2007l ) suggest that the LMC (and with it 
the SMC) are falling into the Milky Way for the first time. 
We find 6 satellites in our ^q" sample that fall into the main 
halo for the first time. Therefore we expect roughly one out 
of the nine dwarfs to be behave likewise. This makes a first 
infall scenario for the LMC not unlikely. 

Leaving aside the three most massive satellites in our 
sample - Sagittarius, the SMC and the LMC - the rest of 
the galaxies share a similar dynamical mass (see Table [T| . 
This makes them useful for probing how environment can 
promote or inhibit star formation. In Figure [TT1 the dwarfs 
are ordered by their star formation histories (c.f. sec. 13. ip . 
We do not see evidence of enhanced star formation at peri- 
cente r, but rather a hint that it is suppressed. This agrees 
with iMaver et al.l (|2007l ). who argue that low pericenters 
are a result of early accretion. Therefore dwarfs that are 
accreted earlier have suppressed star formation because of 



the higher UV background and lower pericenters. However, 
Gaia quality data will be required to test this convincingly. 

Our recovered mean pericentre for all nine dwarfs is 
higher than the mean of both our Zg" and z'Iq samples 
from the VLl simulation. This could hint at satellite de- 
pletion by the galact ic disc as was recently discussed in 
iD'Onghia et al.l (|2009l ). However, proper motion errors will 
bias our pericentre measurement to be systematically high 
(see further discussion below). Our derived apocentre dis- 
tances Ta are lower than the mean of the VLl sub- 
halos (the 50 most massive at redshift z = 0). This dis- 
crepancy cannot be explained by a bias due to the proper 
motion errors as this effect has the wrong sign. Interest- 
ingly, however, the mean of our zf^ sample (the 50 most 
massive before redshift z — 10) agrees well with the ob- 
served mean. This lends further support to the idea that the 
Milky Way's d warfs for med early before reionisation (e.g . 
Efstathiouill992 : Barkan a fc Loeblll999l: iBuUock et al.ll200Cll: 
Benson et al.''2002; Kravtsov et al."20 04 iRicotti fc GnedinI 
2005; Dicmand ct al. 2005; Gncdi n fc Kravtsov! l2006l : 



Moore et al.ll2006l : IMaccio' et al.ll2009D . The mean recovered 
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Figure 12. As in Figure [TT1 but for three different samples taking disc depletion into account. 



Tp/ra ratio is higher in the dwarfs than in the simulation 
samples. Analogously to the pericenter distribution this can 
be explained by the bias from the proper motion errors. 

We investigate the influence of satellite depletion by 
a disc using the three subsets ZxQ{rd ~ lOkpc), z^oifd = 
15kpc) and z^^ird = 20kpc) (50 most massive before z = 
10, without orbits that have pericenters < Vd; 31]). These 
explore difference effective disc sizes. Figure [12] shows the 



same plots as Figure[TT]for the the disc depleted samples. We 
find a trend to higher mean pericentres with larger effective 
disc radii, but also the mean apocenter increases. This is 
due to similar r^/ra ratio in all three samples. Excluding 
Sagittarius, which is in the process of being disrupted by 
the disc, the apocenter distributions in 2:10 (rd — 20kpc) and 
•2^10 (''d = 15kpc) are too high in comparison to the current 
dwarf distribution. This suggests that the total destruction 
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of satellites within > 15kpc is too extreme. With the 
current data it is not possible to discriminate between the 
«io sample without disc depletion and zlo (rd = 10 kpc) with 
a disc with radius ra = 10 kpc. 



4 CONCLUSIONS 

We have evaluated how well we can recover the orbits of 
Milky Way satellites in the light of measurement errors and 
model limitations. To do this, we compared orbits in a high 
resolution cosmological simulation of a Milky Way analogue 
with similar orbits integrated in a fixed background poten- 
tial. We find that: 

(i) With current measurement errors, we can recover the 
last apocentrc Va and pcriccntrc to ^^40% ; 

(ii) With Gaia quality proper motion data, we can recover 
ra and Vp to ~14%, respectively, and the orbital period t 
backwards over two orbits. In this regime, we become lim- 
ited by model systematics rather than measurement error. 
In particular, how well we can approximate the shape of the 
Milky Way potential and how strongly orbits are affected by 
satellite-satellite interactions in infalling loose groups. 

(iii) recovering full 3D orbits - the 3D pericentre rp and 
apocentre ra, remains extremely challenging. This is due to 
a strong dependence on the potential shape that changes 
over time. 

We applied our orbit recovery technique to nine Milky 
Way dwarfs with observed proper motions to determine their 
last pericentre and apocentre distances Vp and Va- We found: 

(i) The mean recovered apocentres are lower than the 
mean of the most massive simulation subhalos at redshift 
z — 0, but consistent with the mean of the most massive 
that form before z — 10. This lends further support to the 
idea that dwarfs formed early before reionisation. 

(ii) With the current data a clear relation between star 
formation history and environment cannot be established. 
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